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ABSTRACT 


Wind energy is currently the fastest-growing energy source in the world. 
However, the inherent characteristic of intermittent energy production, due to 
the stochastic nature of wind, still comprises the main drawback of wind 
power. To avoid such problems, various configurations have been 
reccomended in order to reduce output power variation. The paper 
concentrates on performance benefits of adding energy storage system with 
the wind generator in order to regulate the electric power delivered into the 
power grid. Compared with other means of energy storage, the flywheel 
energy storage system (FESS) is the best choice to solve power quality 
problems. In this paper, a FESS associated to a variable speed wind 
generation (VSWG) is investigated by presenting two control strategies 
applied to the storage system equipped with an induction machine; both 
techniques are studied and developed and consist of a field control (FOC) 
and a Fuzzy Logic Control (FLC). Simulation model is established in 
MATLAB/Simulink and comparative results are then reported. 
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1. INTRODUCTION 

In recent years, wind energy has been the fastest growing and most profitable source of renewable 
energy. Nevertheless, the inherent characteristic of intermittent energy production, due to the highly 
fluctuating and unpredictable character of the wind, still remains the major inconvenient to a wind power. 
Due to the rapid increase in the number of wind farms connected to the network, the variable power produced 
has negative effects on the stability and power quality of the connected electrical equipments [1]. To 
overcome this drawback, various configurations have been reccomended in order to regulate the power flow 
between the wind generator and the power grid. A practical solution consists on introducing an energy 
storage element in connection to a wind power. There are several methods of energy storage that can be 
differentiated into two categories [2]; Long term storage holds energy over a duration ranging from weeks to 
a year such as pumped storage hydropower [3], electrochemical [4], and compress air energy storage [5]. 
Short term storage applies to storage over a duration ranging from several minutes to a few days, such as 
superconducting magnetic energy storage [6], capacitance electric field energy storage [7] and flywheel 
energy storage [8, 9]. The ever-increasing amount of attention on electrochemical energy storage compared 
with other storage systems, but, since rapid response is necessary to compensate power variations in short 
period, they are not appropriate to be associated with WT due to their chemical process. The superconductor 
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coils and the super-capacitors are still under development so they are very expensive [10]. However, 
The inertial storage system by its special features such as low cost, high reliability ,good efficiency, large 
energy storage capacity, rapid response and its long lifetime which is similar to the wind generators; remains 
the most suitable system for storing wind energy. [11-13] prove that the FESS presents an interesting solution 
for adjusting production to consumption. Consequently, when there is an increase in the generated power 
compared to the demanded power, the difference is stocked in the FESS via the electric machine that is used 
as a motor. Inversely, when an imbalance occurs in the power system, the process is reversed and 
the flywheel releases its energy and the machine uses as a generator carrying the network. 

The theory of vector control, applied to the IM of FESS in [14, 15] successfully led to a powerful 
tool for its control. However, experience has highlighted some weaknesses of this method against 
disturbances due to uncertainties of the parameters. It becomes important to use a robust control method, 
insensitive to parameter variations, disturbances and non-linearities. Thus, an intelligent controller such as 
fuzzy controller is needed for improving the speed response; it incorporates human intelligence into 
the process control application which can give better dynamic response of the system [16]. In this paper, a 
low-speed FESS coupled to a VSWG is explored. Figure 1 illustrates a schematic diagram of 
the FESS - VSWG combination under study, where the power generation and storage energy can be linked to 
each other via a DC bus. in such setup, the inertiel storage system guarantees the voltage control of the dc 
link, hence making balance between production and consumption energy. The FESS concerned includes a 
low speed flywheel and squirrel cage induction machine, the latter operates in the area of weakening flux, 
therefore permitting operation at rated power, we present another method control to the FESS, i.e. the Fuzzy 
logic Controller with good comparisons between FOC and FLC are given. The VSWG used is based on a 
DFIG where its stator is directly connected to the grid and its rotor is connected to the grid by power 
converters. To study the power transfer between the wind turbine and the grid, we have used an independent 
power control. So, there were two control blocks (FESS control and DFIG control) in the device under study. 
The first is derived to control the energy storage in the flywheel, the second is derived to control the active 
and reactive power exchanged between the network and DFIG. These two blocks may be separately 
controlled. A MATLAB simulation model has been developed to verify the successful working of 
the considered system, in terms of dynamic response and output power smoothing. 

We have organised the rest of this paper in the following way: Section 2 shows the modelling of 
the FESS (Flywheel, Mechanical shaft and IM modeling). The considered methods of control strategies for 
the FESS (FOC and FLC) are described in more detail in Section 3. Subsequently, Wind turbine modeling 
(turbine, DFIG and Power control of DFIG). Section 5 presents results and analysis. Finally, conclusion is 
given in Section 6. 
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Š 
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Figure 1. Wind system — FESS assembly under study 


2. FLYWHEEL ENERGY STORAGE SYSTEM 
2.1. Flywheel modeling 

The FESS functions as an energy reserve which stocks the energy in kinetic shape in a rotating high 
speed mass that is coupled with the drive shaft of electric machine which supplies an electro-mechanical 
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interface between the wheel and the system. While the charging operation, the machine works as a motor to 
accelerate the rotational movement of the flywheel and thus increase the stored energy, and then the flywheel 
stays in stand-by position. Whereas, on demand, the machine acts as a generator and slows down the wheel 
recuperating the stocked energy in the system [11]. To determine the flywheel inertia, we consider a required 
power during At time. Indeed, to stock the nominal power of the IM P,_1,, during At, the energy AE is then 
necessary such as [12]: 


1 
Ep = 5 Jp OF (1) 
AES Ppt (2) 


Combining (1) and (2), we define the necessary value of the flywheel inertia as [11]: 


2 Pn-imM At 2 Pn-īmM At 
Mea eg ee (3) 
A OF Aimas Amin 
OQemax and Qenin are the maximum and minimum speed limits of the flywheel, respectively. At is 
then the time of storage. This limit should be respected; otherwise we may deteriorate the flywheel energy 
storage process. 


2.2. Mechanical shaft modeling 
The mechanical speed evolution of IM-based FESS can be easily calculated with the dynamic 
formula. The simplified model of this formula is described in [17]: 


ar 


Je —— = Teme f O¢ (4) 


Tem (N-m) : electromagnetic torque, f (N-m-s:rad—1) : viscous friction coefficient. 
2.3. IM modeling 


The IM is selected depending on these benefits in terms of simplicity and robustness of the rotative 
parts. The IM modelled in Park reference frame, is described by [15]: 


X41 —yY.X1ı + wş;.X2 + 7X3 +pxskx, 
T 
X2 — W; X1 — Y.X2 — PXskx3 + =X 
° M 
š3| =| FX 7X3 + (Ws — PXs)X4 fi l (5) 
X4 T- X27 (Ws —PXs) x3 =- 
x p X2. X3 — X1X 
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with x = lisa» lsq Pra» Proh] 
M2 
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3. CONTROL STRATEGY FOR THE FESS 
3.1. Field-oriented control for the FESS 
The reference model in the park, with the rotor flux orientation (@g; = 0, ®ar = Ẹr) has the 


following equations [15]: 
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d. : ; k 1 
—Isq = —Ylsa + Wslsa + = Vv 
dt sd Ylsa s'sq.' T, Pra o.Ls sd 
d. 6 : . 1 
at 1s4 = —Ylsqą — Wslsq — pak Pra + oLs Vsq 
d M. 1 
— — — |] — — 
dt Pra Te sd Ty Pra (6) 


d M. 
qt Pra =0= T, 'sq —(Ws — po) Pra 


d 


d M ‘ 
JuL = JG = PL Pra İsa — Tr — fQ 


The reference value of rotor flux is determined by the flux weakening algorithm. It is 
defined as [15]: 


rn if [Ql < Ofn 
Pr-ref Nr) = On (7) 


Orn: if [Ql > Ofn 
|Q¢| 


The FESS reference power P+ ref should be restricted to the nominal value of IM power to avoid the 
overheating IM. The reference torque is indicated as follows: 


Pt_ref 


Tim-ref = Qf (8) 


The quadratic reference current becomes: 


TīM-ref -Lr—IM 
P.M.Ọrd-ref 


(9) 


İsq-ref = 
The reference power exchanged between the storage system and the DC bus is given as follows: 
Pi-ref = Preg z Pwg + AP (10) 


Where, Preg the reference grid active power, and Py, the power generated by wind generator, AP the 
power fluctuation on bus capacitor [17]. 


and, the kinetic energy E,;¢¢ is computed by the integration of (10) [11]: 
t 1 
Ecref = iF Pr ref dt + Eco = J Je OF (11) 


where Eo the initially kinetic of the storage system [11]. The flywheel speed reference is expressed as: 


2 Ecre 
Ores = J (12) 


3.2. Fuzzy logic control for the FESS 

Figure 2 illustrates the block diagram of fuzzy controller, includes four parts: - Fuzzification - 
knowledge rule base - fuzzy inference - Defuzzification [16-18]. 

The conception of the fuzzy logic controller begins by allocating the input and output variables. 


The more important variables entering the fuzzy logic speed controller where selected as the speed error and 
its temporal variation. Two inputs variables EQ(k) and dEQ(k), are computed at each sampling time [19]: 


EQ(k) = Afrek) — Of Ck) (13) 
dEQ(k) = EQ(k) — EQ(k—-1) (14) 
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Where E 0(k — 1) is the error value at the preceding sampling instant. The output variable of fuzzy 
logic speed controller is the varying control current disqe (K) that is incorporated to obtain the reference 


control current, i (k), as indicated in the next formula. 


Sdref 


İsqrer K) - İsqrer K 7 1) T disqrer (K) (15) 
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Figure 2. Fuzzy Control system scheme 


3.2.1.Fuzzification 

The quality of this step is critical to the success of this work. At this phase, the crisp variables 
EQ(k) and dEQ(k) are respectively transformed in to fuzzy variables EQ and dEQ. Triangular fuzzy 
membership functions were defined as shown in Figure 3. The discourse universe of all input and output 
variables is set as (-0.8, 0.8). Appropriate scaling coefficients are selected to bring the input and output 
variables into this universe of discourse. Each discourse universe is split into seven overlapping fuzzy 
clusters: NL (Negative Large), NM (Negative Medium), NS (Negative Small), ZE (Zero), PS (Positive 
Small), PM (Positive Medium) and PL (Positive Large). Each fuzzy variable is a member of the subsets with 
a degree of membership u ranging from 0 (non-member) to 1 (full member). All membership functions have 
an asymmetric shape with more clutter near the origin (stationary state). This allows a greater precision in the 
steady state [19]. 


(EQ), (AEQ), u(Adisq,.) i E NL NM NS ZE PS PM PL 











NL NM NS ZEPS PM 


EQAE, Adige PL ZE PS _ PM_PL PL PL PL 


Figure 3. Membership inputs/outputs 


3.2.2.Knowledge base and Inference engine 

The fuzzy rules are expressed using the IF-THEN form as follows: if (input! and input2) then 
output, in example: IF (EQ is Positive Small) AND (dEQ is Negative Medium) THEN (digg, is Negative 
Small). TThe fuzzy variables EQ and dEQ are processed by an inference mechanism based on a set of control 
rules contained in (7*7) table as shown in Table |. The crisp output of the FLC is obtained by using MAX- 
MIN inference algorithm [19]. 


3.2.3.Defuzzification 

In this step, a sharp value of the output variable di,, -(k) is determined through the defuzzification 
method, which first evaluates the centroid of every output membership feature for each rule. Then the final 
output is computed as the mean of the single centroids, which is weighted by their heights (degree of 


membership) as shown below: 
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16 
Ts (disaren) a 


disg,o¢(K) — 


The control system scheme of the system is detailed in Figure 4. 
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Figure 4. Fuzzy Control system of inertial storage scheme 
4. WIND TURBINE MODEL 
4.1. Turbine modeling 
The modeling of a wind energy system is given by [20]: 
3 
P, _ P ~~ 
Cy(A,B)p S v3 Paar Cp(A,B)p S v3 (17) 
Paer = De and Taer = “OL = ar a 


The Power coefficient C, dependent on the blade tilt angle B and the tip speed ratio (TSR) A has the 
following definition [20, 21]: 





21 
C A, B) = 0.5176 a —0.4B-5) e™ + 0.0068 À 
1 = 1 0.035 gy a HR (18) 
i = ooog Bi POLT | 


Where P,, Pàers Taer, P, S, R, Qi and v are the turbine power, the aero-dynamic torque, the aero- 
density (p = 1.22 kg / m3), the circular area swept by the turbine, the turbine radius, the turbine speed (rad / s) 
and the wind speed (m/s), respectively. The turbine is connected via a gearbox to the generator; the gearbox 
is modelised with these two equations [22]. Where: Tmec 1s mechanical torque, Qmec 18 the generator's speed, 
G is the ratio's gearbox. 





T ur 
Tinet = F 
Omec = GO (19) 
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The basic dynamic equation of mechanical system on a shaft of the DFIG is as follows [20-22]: 
where: 
f : Coefficient of friction, Tem : electro-magnetic torque. 


dOmec 
a Fi Unec Imee Tem (20) 


4.2. Model of doubly fed induction generator (DFIG) 
The Park state model of the generator is represented in the following matrix form [23, 24]: 


: d . d 
Vsa = Rglisa + at Psd — DsPsq Vra = Rrira + at Prad — Or Org 
(21) 


; d : d 
Vsq = Rsisq + at Psa + WsPsq (Vig = Relrg + qe Pra PrPrd 


Rs R; : stator and rotor phase resistances. ws = œ + wr, and w = P. Qmec is the electrical speed P 

is the pair pole number. The stator and rotor flux takes the following form: 
vi = Ls-İsa + M. ira va = Leira + M [isa (22) 

Psa = Ls-isq +M .irqa (Pra = Lr-irg + M. isg 


In which, iggy sg /İra»İrq are d.q components of the stator and rotor currents, Ls, L}, M: are stator, 


rotor and mutual inductances. The stator active and reactive power and the electro- magnetic torque of the 
DFIG takes the following form [23]: 


A = Vzd- isd + Vsq- İsq (23) 
Qs = Vsq: İsa — Vsa- İsq 

3 M ; : 
Tem = > P Ls (@sq- ira E Psa- İrq) (24) 


To facilitate the control of electricity generation from wind, an independent control of active and 
reactive power will be achieved by developing the equations that relate the values of the rotor voltages 
generated by an inverter to the stator active and reactive power. [25]. Assuming a constant flux in the stator, 
we obtain: 


iy = 0 (25) 
Vas = Vs 

to one (26) 
Qs = Vs las 


The relationship between the rotoric and statoric currents is: 


oy ke. M. Ps 

fas F 
M. (27) 
lqs = — ee tar 


To obtain the expression of powers in function of stator currents, we substitute the currents in (26) 
by (27): 








= M : 
Ps = — Vs = lar 28 
Q, = Veo =. VsM į (28) 
$2 Ls@s5 Ls dr 


In order to correctly controlling the machine, we need to find the relationship between the rotor 
current and voltages: [26] 
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Var = Relar + (br E)E -go (Lr E)r 





(29) 








Where g is the slip of the asynchronous machine. A complete DFIG model scheme is 
depicted in Figure 5. 


4.3. Indirect field oriented control (IFOC) of the DFIG 

To ensure stable functioning and to allow independent control of the active and reactive power of 
DFIG, a template-based PI regulator is constructed utilizing the dynamic equations. A schematic diagram is 
given in Figure 5 [25]. 
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Figure 5. Control scheme of DFIG 
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5. RESULTS AND ANALYSIS 

The parameters used can be found in Tables 2 and 3. The simulations were performed in Matlab / 
Simulink environment. Our aim is to show the behavior of the FESS on a charge/discharge cycles using 
speed control strategy. Figure 6 shows simulation results achieved by analysing the efficiency of FESS with 
IM. In this case, the FESS is not linked to the VSWG. Then, the flywheel decelerates and discharges its 
energy into the grid via the IM. At speed greater than 157 rad/s, the IM connected to the flywheel functions 
in the field weakening zone. For 0 <t< 1.5 s and P;-ref > 0, the flywheel is accelerated to its maximum speed 
and stocks energy in the grid. It functions in motor operation. However, for 1.5 <t <3 s and Pr_rer < O, the 
flywheel decelerates and runs in generator mode. To test the robustness of the parametric fuzzy controller, we 
performed a parametric variation of the machine and studied the influence on the performance of the speed 
control and its power comportement. The variation is plotted on the rotor resistance Rr (50% increase of Rr). 
The results achieved Figure 7, indicate that the variation of the rotor resistance influences the decoupling of 
the system for PI regulator; however, the fuzzy control almost keeps the original performance despite the 
variation of parameter in question. The simulation of the indirect field oriented control of a DFIG-based wind 
system is represented in Figure 8. This figure demonstrates that our system has satisfactory dynamic, the 
active and reactive powers of the stator follow their references almost perfectly and provide a perfect 
decoupling between both axes. 


Table 2. FESS Parameters Table 3. DFIG and Turbine Parameters 
P=2 R, = 4.850 R, = 3.8050 P, = 10KW R, = 0.4550 G=90 
L, = 0.274 H L- = 0.27 Hen M = 0.25 Hen R, = 0.1900 L, = 0.070H R= 35m 
J; = 20.28 m2kg fe = 0.022 Nm/s L, = 0.0213H M=0.034H N=3 


In what follows, we passed to FESS - VSWG connection in order to regulate the electrical power 
generated by the wind. All simulations are run with the same wind profile as shown in Figure 9. The rated 
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power of the DFIG is 10KW and the IM was rated at 3 kW at 1500 Tr/mn, the variable power of the wind 
turbine and the desired grid power (fixed to -2.2 kw) are shown in Figure 10. Figure 11 shows the power 
transmission between the FESS and the grid, i.e power storage of the inertial system (FESS). The flow 
direction of this power fluctuates as a function of the power generated by the wind turbine. If the latter is less 
than -2.2 KW, the power transfer from the network to the SISE (storage); if not, from the SISE to the 
network (destocking) to offset for the power deficit of the wind. A rotational speed increases when the energy 
is transferred to the FESS, and decreases when the energy is transferred to the grid. Figure 12 shows the 
electro-magnetic torque of the IM which is fluctuant; it depends directly on the electrical power of the FESS, 
which is reflected in its shape identical to that of the power of the FESS. The direct and quadrature rotor flux 
elements of the IM are also depicted in Figure 12. The quadrature component is always zero, justifying the 
rotor flux-oriented control. 
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Figure 12. Electromagnetic torque and rotor flux of the IM 


The comparison between both controls reveals several benefits of the fuzzy logic controller, such as 
robustness and better reference tracking compared to the flux oriented control. In the paper [14], it is 
apparent that the FESS power exhibits frequent oscillations, implying a lower robustness. When comparing 
these results with ours, we find that the Flywheel energy storage system in our study is more reliably and 
robust, which allows for better energy injection into the network. 


6. CONCLUSION 

The research reported in this paper is focused on the analysis, modelling and simulation of a 
variable speed wind turbine using a double-fed induction machine in combination with a flywheel storage 
system in the aim to resolve the problem of fluctuating power output. Two FESS supervision strategies 
namely FOC and FLC have been applied. It has been noted that the FLC is a preferable choice for such 
application. The simulation results displayed prove that a correctly functioning of the storage system is in 
place, in effect, the measured power exactly tracks the reference power. The good control effect is proved by 
the simulation results. The process of charging and discharging are stable and the control 
objectives are achieved. 
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